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Abstract 

We show that critical opalescence, a clear signature of second-order phase transition 
in conventional matter, manifests itself as critical intermittency in QCD matter pro- 
duced in experiments with nuclei. This behaviour is revealed in transverse momentum 
spectra as a pattern of power laws in factorial moments, to all orders, associated with 
baryon production. This phenomenon together with a similar effect in the isoscalar 
sector of pions (sigma mode) provide us with a set of observables associated with the 
search for the QCD critical point in experiments with nuclei at high energies. 
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It has been suggested recently that the baryon number fluctuations, in event-by-event 
studies of multiparticle production, may lead to new observables, in experiments with nuclei, 
associated with the existence and location of the critical point in the QCD phase diagram 
[1]. Theoretically, the significance of the baryon number density ^^(x) near the critical point 
comes from the fact that its fluctuations, expressed in terms of the density-density correlator 
(n_B(x)nB(0)), follow the same power law as the sigma fleld correlator {a{x)a{0)) associated 
with the chiral order parameter (t{x) ~ [1)2]- As a result, the critical thermodynamics of 
QCD matter in equilibrium can be formulated in terms of the baryon number density ^^(x), 
an alternative but equivalent order parameter related to the final state of nuclear collisions 
when they freeze out close to the critical point. The corresponding theory of critical power 
laws in the baryon sector, the observability of which we explore in this Letter, is based on 
the effective action belonging to 3d Ising universality class: 
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where 5 ~ 5 is the isotherm critical exponent and G a universal, dimensionless coupling in 
the effective potential, with a value in the range G ~ 1.5 — 2 [3]. In writing eq. (1) we 
have considered, following the above discussion, the order parameter m{x) — gT~'^nB (x), 
associated with the fiuctuations of baryon number density at the critical point; the factor 
T~'^ guarantees dimensional consistency and g is nonuniversal dimensionless constant. 

The free energy (1) in order to describe a critical system of QCD matter produced in 
nuclear collisions, must be adapted to the rclativistic geometry of the collision. To this end, 
the longitudinal coordinate a:|| is replaced by the space-time rapidity y and the corresponding 
measure of integration in (1) takes the form: dx\\ — TcCOshydy where Tc is the formation 
time of the critical point. In the central region of size Sy the configurations ub {y, x±) 
contributing to the partition function are boost invariant quantities defining at the same 
time two dimensional baryon-number density configurations Pb {x±) in the transverse plane: 
n'B{x±) = Pb {x±) [2rcSmh{6y/2)]~^ . Integrating now in rapidity and rescaling the basic 
variables: x± — > TcX±, pb T~'^pB, eq. (1) is simplified as follows: 



r, Ipb] = Gg' J d'x^ I VxPbI' + G {gG)'-' p%+' 
G = {%Teff)~^ , Te// = 2rc siuh {5y/2) 



(2) 
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The critical fluctuations of systems belonging to the class (2) of a 2d effective action can be 
consistently described in a scheme where the saturation of the partition function is obtained 
considering the contribution of the singular solutions p^^ of the Euler-Lagrange equation [4] : 



Z ^ ^ e-^^I"^^'] ; V'J$ -{5 + l)G {gCf-' [p^^^\ " = (3) 

s 

The description of the critical system in this scheme is optimal if Cg'^ 1. 

The main characteristics of critical QCD matter (5 ~ 5, G ~ 2) produced in nuclear 
collisions and described by eqs. (3) are summarized as follows [4]: 

(a) The system is organized within critical domains (clusters) of a maximal size (correlation 
length): ^± ~ ^ sinh ((5y/2). 

(b) For typical nuclear collisions (tc ~ 10 fm, Sy > 2) the transverse correlation length 
becomes sufficiently large > 6 fm) allowing for critical fluctuations to develop in full 
strength when the system freezes out near the critical point. 

(c) Within a critical domain (|(^x_l| < ^_l) the correlator obeys a power law corresponding 
to a fractal dimension dp = {dp ~ |): (pb(^±)pb(0)) ~ l^j,!"^^"^. This power law is the 
origin of critical opalescence [5] in baryonic QCD matter produced in high-energy nuclear 
collisions. 

In conventional matter (QED matter) the phenomenon of critical opalescence is revealed 
when light of long wavelength (comparable to the correlation length) scatters from a sub- 
stance near criticahty. The intensity of scattered hght is proportional to the Fourier trans- 
form of the density-density correlator and becomes singular for small momentum transfer k 
giving rise to a macroscopic effect [5] . In QCD matter, correspondingly, the Fourier transform 
of the correlator {pb{x±)pb{0)) obeys also a power law: 

PkPl) - J d'x^e-^'-'- {pB{x^)PBm 

{pB{x^)pBm ~ \x^f--' ; ^ (4) 
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where is the Fourier transform of pb{x±) and \^PkP\^ the two-particle correlator in trans- 
verse momentum plane, associated with baryon production in nuclear collisions. More pre- 
cisely, 

[pkPu) = E (PBip)PBip + k)) ~ {pB{0)pB{k)) (5) 



under the assumption that the dependence oi (^pb{p)Pb{p + k)j on the reference momentum 
p, is weak. 

The power law (4) in momentum space can be generahzed to a self-similarity property 
of multiparticle correlators, valid for the solution (3) of QCD matter at criticahty [4]: 

PB{Xh) ■ ■■Pb{XK-i)Pb{0)) = A-'^-(''-^) (pBih) ■ ■ ■Pb{K-i)Pb{0)) , (g = 2, 3, . . .) (6) 



In geometrical language, eqs. (6) define a fractal dimension dp in transverse momentum 
space, independent of q {dp — 2 — dp), directly related to the isotherm critical exponent of 
QCD: dp — [6]. The observable effect imphed by (6) is a specific class of power laws, 
satisfied by the scaled factorial moments of all orders [7] : 

-{NB{NB-l)---{NB-q+l)y 



F,{M) = 



{q — 2,3...) defined in small domains 5Sm of transverse momentum plane, constructed by 
a subdivision of the available space in two-dimensional cells {6Sm ~ M~^). The power 
laws (7) with a characteristic linear spectrum of indices = (2 — dp){q — 1), describe, in 
general, the fluctuations in a second-order phase transition [6,8] but here, in particular, they 
reveal the effect of critical intermittency associated with the production of baryonic QCD 
matter near the critical point. This phenomenon is the analogue of critical opalescence, 
observed in conventional matter near criticality. Both phenomena have their origin in the 
power laws of the correlators in configuration space, implying giant density fluctuations at 
the critical point as a result of the appropriate fractal geometry, valid within the universality 
class of the critical system under consideration. 

In order to make the prediction (7) a transparent signature for the QCD critical point in 
experiments with nuclei, one may proceed to a Monte Carlo simulation of critical events in 
the case of a typical process A-l-A at high energies. We have chosen collisions of medium size 
nuclei producing net baryons in the central region with average multiplicity {Nb) ~ 60, at 
energies corresponding to CERN-SPS. At the critical point we expect an extended plateau 
in the rapidity spectrum of net baryons [1] allowing for a choice (5y ~ 5 in the system under 
consideration. A self-consistency requirement in this treatment is the constraint Cg'^ ^ 1 
which may justify a posteriori the saddle point saturation of the partition function leading 
to eqs. (3). In order to clarify this issue for the above system under simulation, we consider 
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the following generic form coming from the thermal average {ps) with the aid of eqs. (3): 

(^^5, G^2, r,^10/m) (8) 

where Scr is the area of the critical cluster, Scr — For the specifications of the system 

under study we find Cg^ ~ O(IO^), a result which fulfils comfortably the consistency 
requirement {Cg^ ^ 1). 

The simulation of events involving critical fluctuations in the baryonic density requires 
the generation of baryon transverse momenta correlated according to the power-law (4). 
One possibility to produce such a set of momenta is to use the method of Levy random 
walks [9] . It is convenient to perform independent one-dimensional walks in each transverse 
momentum component separately and then form the Cartesian product in order to obtain 
the corresponding vectors. In this case the successive steps in each dimension are chosen 
according to the probability density: 

p{Pi) = N'^ A"^"" i Pmin < Pi < Pmax, i = X,y (9) 

■j^ I Pmin \ 

\Pmax / 

To generate a single event one performs n — 1 random walk steps in each direction {p^ or 
Py) where n is the multiplicity of the event. At each step the updated position of the walker 
determines the Px (or py) coordinate of the baryon transverse momentum respectively. The 
starting transverse momentum vector in each event is chosen uniformly in [—Pmax,Pmax] x 
[— Pmax, Pmax]- The produccd set of baryon transverse momenta possesses the correct fractal 
dimension dp = \ provided that ^ 10^ and u = ^. The described algorithm can be 
used to generate an ensemble of events involving critical baryon density fluctuations. The 
corresponding fractal pattern can be revealed through factorial moment analysis [7]. In 
Fig. 1 we show, in a log-log plot, the results obtained for an ensemble of 100 critical events 
each one having the multiplicity n = {Nb) = 60. The full circles are the calculated moments 
up to the 4th order. The solid lines arc the theoretical predictions according to eq. (7). It 
is clearly seen that the factorial moments follow to a good approximation and for a wide 
range of scales a power-law behaviour Fg ~ (M^)*« with exponents Sg satisfying very well 
the theoretical prediction Sq — {q — — ^) . In the same plot we show for comparison the 
behaviour of the second moment (grey triangles) in a conventional system corresponding to 
the formation of mixed events from the original ensemble of the critical events. As expected. 
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the effect of critical fluctuations has disappeared in the system of mixed events (the slope 
st^ ~ 0). 




Fig. 1. The log-log plot of the factorial moments of order q = 2,3,4 (full circles) for 
an ensemble of 100 critical events generated by the Monte- Carlo algorithm described in the 
text. The corresponding theoretical power-law predictions are shown with solid lines. The 
grey triangles present the result for the second moment for a set of mixed events. 

We have shown that baryon-number density fluctuations, near the QCD critical point, 
develop a distinct pattern of power laws, in transverse momentum plane, associated with the 
effect of critical intermittency in QCD matter. The observability of this effect is enhanced 
by the conjecture that the same power laws (Fig. 1) arc expected also in the case of net 
proton-number density fluctuations, avoiding therefore the observational ambiguities from 
the contribution of neutrons to the factorial moments (7). In fact, it has been argued by 
Hatta and Stephanov [1] that the singular parts of baryon susceptibility {xb) and proton 
number fluctuations coincide, xb ~ (^^p-p^^p-p) • As a result, the power laws associated 
with these singularities are the same for both the baryon-number and proton-number density 
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correlators. With these remarks we may conclude that proton number measurements in 
transverse momentum plane may reveal the effect of critical baryon-number intermittency in 
nuclear collisions, as a signal of QCD criticality. If one combines this observable effect with 
the corresponding phenomenon in the sigma mode of pion production [10] then a complete 
set of observables, associated with the existence and location of the QCD critical point, 
naturally emerges. 
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